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Abstract. Given a non-oscillating gradient trajectory I7I of a real an- 
alytic function /, we show that the limit v of the secants at the limit 
point O of | 7 | along the trajectory \^\ is an eigen-vector of the limit of 
the direction of the Hessian matrix Hess(/) at O along |7|. The same 
holds true at infinity if the function is globally subanalytic. We also 
deduce some interesting estimates along the trajectory. Away from the 
ends of the ambient space, this property is of metric nature and still 
holds in a general Riemannian analytic setting. 



1. INTRODUCTION 

Let / : (M n ,0) — > (R, 0) be the germ at O of a real analytic function 
singular at O. The behavior of the germ at O of gradient trajectories of 
the function / with cu-limit point O is barely understood. Thorn's Gradient 
conjecture, proved by Kurdyka, Mostowski and Parusihski [S], states that at 
O the trajectory has a tangent. In whole generality nothing else is known. 
The plane smooth case is very well understood, and (bounded) plane gra- 
dients curves are definable in the Pfaffian closure of M an . But a few special 
cases trying to address the non-spiraling/non-oscillating behavior of such 
trajectories [1 3|. HI [6], and the main result of [12], that is all that is known. 

The aim of this note is to report an intriguing and not expected result, 
along gradient trajectories that are non-oscillating at their limit point. We 
recall that a gradient trajectory 7 is (semi- analytically) non- oscillating at 
its limit point O € crit(/) if, given any semi-analytic set H (germ at 0), 
either j<zHorjnH = @ (as germ at O). Our main result is 

Theorem 15.51 Let 7 be a gradient trajectory of a real analytic function 
f : (W l ,0) — > (R, 0). Assume that 7 is non-oscillating at its uj -limit point 
O € crit(/). Let v € S n_1 be the limit at O of the secant along 7. 

1) The oriented direction of the matrices Hess(/) has a (non-zero) limit Ti 
at O along 7. 

. , Hess(/)(x) • V/(x) 

2) hm — 7—7-7 — z — =-rz , = —v. Thus the unit vector v is an 

J | 7 | 9x _>0 |Hess(/)(x) • Vf(x)\ 

eig en- direction ofH. 
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This result is of metric nature, and although it is here just stated for the Eu- 
clidean metric it still holds in any analytic Riemannian setting (see Section 
[8]). This result is somehow providing a clue along non-oscillating gradient 
trajectories about where a gradient trajectory is ending-up. Note that this 
result was not even known in the plane case, where it applies. 

The paper is organized as follows. 

Section [2] introduces the few useful definitions, notions and notations we 
will deal with. 

In Section [3] we present quickly the framework which we will operate in, 
and recall [§] an interesting technical result for us. 

Section [J] is aimed at giving motivations to this paper, in stating and 
re-proving results about links between gradient trajectories and Hessian 
matrices that should be well known folklore. These special cases suggest 
two possible different interpretations of these connections between Hessian 
matrices and gradient trajectories. 

The main result is stated and proved in Section [5j 

It is followed by a few corollaries in Section [6] presenting some estimates 
along the trajectory. 

Section [7] deals with the analog result at the limit point at infinity of a 
non-oscillating (at infinity) unbounded gradient trajectory for of globally 
subanalytic real analytic function. We also provide the corresponding esti- 
mates expected from those found in the previous section. 

Section [8] we state and prove our main result in the general (real analytic) 
Riemannian case. 

Thanks. I am very grateful to the Fields Institute for its support and 
the working conditions provided since this work started during the Fields 
Institute's Thematic Program on O-minimal Structures and Real Analytic 
Geometry January-June 2009. I thank F. Sanz a lot for having suggested 
to jump from the plane case to the non-oscillating case. I owe many thanks 
to K. Kurdyka for precious comments, suggestions and encouragements. 

2. Background - Notations 

Let 7 : 1 — > W l be C 1 parameterized curve over a connected interval I with 
non-empty interior. We will write | 1 for the support of the parameterized 
curve, that is 

\M = {<t),tel}. 

Definition 2.1. The curve is non-oscillating if for any semi-analytic 
subset H, the intersection H n I7I has finitely many connected components. 

Since we are interested in a local problem we prefer using the following 
local version: 

Definition 2.2. Assume that the origin O belongs to clos ( 1 7 1 ) \ |7|. The 
curve I7I is non- oscillating at O, if for any semi-analytic subset H , the 
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germ at O of the intersection H n I7I is either empty or is the germ of | — y | 
at O. 

These notions were explored for real analytic vector fields |131 [21 [1] on 
3-manifolds to study their asymptotic dynamics at their singular points in 
certain situations, namely to discuss spiraling. 

The pointed vector space Sym(IR n )* of non-zero real symmetric matri- 
ces of 1™ is an open cone on the unit sphere of Sym(R n ) with vertex the 
null matrix, on which the group (M>o, •) acts naturally, smoothly and semi- 
algebraically by the homotheties of positive ratio. The resulting quotient 
space Sym(n) is the space of oriented directions of the real symmetric n- 

matrices and is a real algebraic manifold isomorphic to S 2 _i . The 
quotient mapping a : Sym(R n )* — > Sym(n) is a smooth semi-algebraic sub- 
mersion. 

Let / be a C 2 function defined on some open set U of W 1 equipped with 
some coordinates. We define, Hess(/) = [Kj] to be the Hessian matrix of /, 
and the function h on U as h{x) = j hij( x ) 2 ] 1 ^ 2 - We wn ^ write H(f) for 
the matrix /i _1 Hess(/) where it is defined. We obviously get <r(Hess(/)) = 
a(H(f)). Note that H(f) is not a(H(f)), but has a similar property: let 
(Hess(/)(xfc))fc € Sym(R n )* be a sequence converging to the null matrix. 
The accumulation values of the sequence (H(f)(xk))k are bounded and will 
not all be the null matrix. In particular if the sequence (H(f)(xk))k has a 
unique accumulation value T~L, the sequence cr(H(f)(xk))k then converges to 
cr(7i) = XH for some positive real number A. With an obvious but small 
abuse of language we will also call such a limit % an oriented limit direction 
of the sequence (H(f)(xk))k- 

Last a short word on notations. 

Unless explicitly mentioned M. n comes equipped with the Euclidean struc- 
ture and we will write {—,—) for the scalar product and | — | for the associated 
norm. 

3. Setting 

Let / : U — > K be a real analytic function defined on an open subset 
U C M. n containing the origin O, such that O S crit(/) the critical locus of 
/. We also assume f(0) = 0. 

Let ut be the unitary gradient vector field defined on R n \crit(/). But its 
orientation, this vector field does depend only on the foliation by the con- 
nected components of the levels of the function, not on the chosen function. 

Let us consider the following gradient differential equation: 

j(s) = 17(7(3)) with 7(0) = x and f(x ) < 0. 

Given any maximal non stationary solution 7, let I7I be the underlying 
subset of the parameterized integral curve, namely 
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h\ = Ms),s>o}. 

Definition 3.1. We call \ j\ a half-trajectory of the gradient field associated 
with the function f. 

For short we may say trajectory instead of half-trajectory. 

Assume we are given a half-trajectory I7I, we assume that the function 
s — > fi'y(s)) increases to and that the omega limit set w(|7|) contains O. 
Let l(\j\) be the length of [7] between 7(0) and the origin O. We know since 
Lojasiewicz that it is finite [11], and so w(|7|) = {O}. 

The next result is the first systematic result on the behavior of gradient 
trajectories at their limit point, and so far the only such one within the 
category of analytic functions. 

Theorem 3.2 (|8j). The length of the radial projection t — > onto S n_1 
of the gradient curve 7 is finite. 

The straightforward and looked for consequence of Kurdyka, Mostowski & 
Parusihski's result is that the well-known Thorn's Gradient Conjecture holds 
true: 

Thorn's Gradient Conjecture. There exists v E S n_1 such that 

Tt S 1 X 

lim - — — - = v or equivalently lim - — - = v. 

\j\3x-+0 \x\ 

Writing d r for the radial vector field, the gradient vector field V/ then 
decomposes into two orthogonal components as V/ = d r fd r + V/ where 
d r fd r is its radial component and V/ is its spherical part. 

To prove their result, Kurdyka & Als used fine estimates along the tra- 
jectory. In particular we mention the following result that will prove to be 
useful for our purpose: 

Proposition 3.3 (|8j). 1) There exists a rational number m := m(|7|) and 
a positive real number a := a(|7|) such that along 

f = — ar m + o(r m ) and d r f = —mar" 1 ^ 1 + o(r m " 1 ) 

2) There exists a finite subset S C Q>o x IK>o of pairs (l,b), such that for 
any trajectory a, such that w(|a|) = {0} and f is negative along \a\, the 
pair (m(\a\) , a(\a\) as in point 1) corresponding to \a\ belongs to S. 

4. Hessian and Gradient: special cases 

This section is devoted to present some examples exhibiting strong con- 
nections between limits of direction of Hessian matrices and gradient vector 
fields. This material should be well known and at least part of the folklore 
of this subject. There are somehow similar suggestions in |12| . 

1) We decompose the real analytic function / as the sum of its homoge- 
neous parts at O, namely 
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/ = f p + f p+1 + . . . 

Let Hess(/) be the Hessian matrix of /. 

We recall that for each k > p the following formulae hold 

k(k - l)f k = (k- l)(Vf k ,rd r ) = (Hess(A) • rd r ,rd r ), 

Let j3 : S n_1 xlR + — > W 1 be the spherical blowing-up with center O defined as 
(u, r) — > ru. Up to a multiplication by p(p — l)r 2 ~ p the gradient differential 
equation now reads: 

r = r (Hess( / p ) (u) ■ u,u) + 0(r 2 ) 

u = Hess(/ P )(u) • u — (Hess(/ P )(n) • n, u)u + 0(r). 

The vector field corresponding to this differential equation is called the di- 
vided gradient vector field, and is denoted £/. It extends analytically on 
S"" 1 x R. 

Thus if v € S n_1 is not an eigen-vector of Hess(/ P )(i;), then u\ r= o jU=v ^ 0. 
Thus any limit point v of a gradient trajectory is such that v is an eigen- 
vector of Hess(/ P )(^). (Which is not such an interesting information once 
Hess(/ p )(^) is the null-matrix.) 

The following elementary result should be known, although we never came 
across it this way. 

Proposition 4.1. Let f be a homogeneous polynomial of degree p > 2. 
Then along any gradient trajectory \j\ with oj-limit set O with limit of secants 
v, we find f(v) ^ so that v is an eigen-vector of the non null-matrix 
Hess(/)(^) associated with a non-zero eigen-value, and so the vector field 
is reduced at v, that is has a non-nilpotent linear part. 

Proof. Let I7I be a trajectory of V/ with w-limit point O and limit of secants 
v. Let r(s) = |7(s)| and let u(s) = ^(7(5)), then / 07(5) = r(s) p f(v(s)). 
Proposition 13.31 gives along I7I, rd r f = mf + o(r m ) = —mar" 1 + o(r m ) as 
r(s) —> 0. Since rd r f = (Vf,r5 r ) = pf, we must have p = m and so 
/(*/)- -a /()). 

The statement on £y holds since (Hess(/)(z^) • v, v) = p(p — l)f(v). □ 

The following result should also be already known. Nevertheless, we state 
it and provide a proof. 

Proposition 4.2. Let f = f p + f P +i + ■ ■ be the sum of the homogeneous 
parts of the real analytic function f . Let v be the limit of the secant of a 
trajectory (7! of the gradient V/ with w(|7|) = {O}, and such that f p {v) = 
—a ^ 0. Then the divided gradient vector field £j is resolved at u, that 
is has a non-nilpotent linear part at v. Moreover, the eigen-values of 
the linear part of are entirely determined by Hess(/ P )(^) (which is also 
^x/\x\->u [x|- p+2 Hess(/)(x)j. 
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Proof. Let us look at the dominant part of the gradient vector field nearby 
the oriented direction v. First we assume that v = d Xl . Second if H := 
Hess(/ P )(^), and since v is an eigen- vector of H, we take orthogonal coor- 
dinates x = (xi,y) at so that H = [hi] is diagonal, where hi is the i-th 
diagonal coefficient and hi = —p(p — l)a ^ 0. 

Since u\ = J 1 — Yl v i > where v i = Ui for i > 2, we find that 
(Hess(/ P )(u) -u,u) = ht + ^J(^i - h)vf. 

i>2 

Thus, in a neighborhood of v the divided gradient differential equation be- 
comes 

, , I r = h ir + 0(\(v,r)\i) 

1 J \vi = - /ii)^ + 0(H 3 ) + 0(r), 

As announced we find that the eigen- values of the linear part of £y are easily 
deduced from that of Hess(/ P )(i^). □ 

The following is a straightforward consequence of Proposition 14.21 

Corollary 4.3. Under the hypothesis of Proposition \4^M we find 

lim H(f)(x) = CRess(f p )(u), 

|7|3:r->0 

for a positive constant C . 

The second point we would like to put in light is the following: 

2-1) Assume that / = f p is a homogeneous polynomial of degree p > 3, let 
£ be the set of points where the gradient V/ is an eigen- vector of Hess(/), 
that is 

E = {x : df A d|V/| = 0} C {x : d/ A d|V/| 2 = 0}. 
It is a closed semi-algebraic cone over O. 

In particular each open ray M>oit, for u £ (E n S n_1 ) with /(«) < 0, is a 
gradient trajectory with w-limit point O. 

For a positive number r, we define E r = E n f~ 1 (r). It has finitely many 
connected components on each of which |V/| takes a constant 

values pcj(r). Note that U r >oE* is a positive open cone C* := M>oE* , 
where E l is a closed semi-algebraic subset of S ra_1 . At any x G C*, using 
the spherical coordinates and writing x = ru, we get 

|V/(tt)| = pr x ^ p Ci{r) = pa. 

Up to a re-numeration we have < c\ < . . . < c q . At any x = ru G C*, we 
have Hess(/)(x) • Vf(x) = \{x)v{x). We deduce that A(x) = r p_2 A(u), and 
so the mapping T, 1 3 u ^ X(u) is semi-algebraic. 

If v is a limit of secants at O of a gradient trajectory (7!, then we know 
that Hess(/)(^) • u = \v for A 7^ 0. So \j\ = M>qi / and is contained in E. 
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Since | — y | C C* , we get that V/ o j(r) = —pr p 1 civ, since we have agreed 
that / o 7 < 0. So we find that 

\(u) = p (p - l) Ci and so f{v) = -Ci. 

This very special example suggests the existence of links between gradient 
trajectories, Hessian matrices and their eigen- values and limits of extrema 
of |V/| restricted to the foliation by the level of /. 

2-2) Let / be a real analytic Morse function, vanishing at the origin. We 
write / = /a + 0(|x| 3 ). The divided gradient differential equation ensures 
that the limit v of the secants of any trajectory \^\ falling onto O is an 
eigen-direction of the Hessian matrix Hess(/2), which is invertible at O. Up 
to a an orthonormal change of coordinates we write 

/(*) = - £-=i ^ + EU n*i + °(M 3 )> 

where s, u are non negative integers and s + u = n and, when s > then 
< Ai < . . . < X s and, when u > then < fj,± < . . . < fi u . From Corollary 
14.31 along a trajectory | — y | we get r _1 7 — > v and that v lies in the eigen- 
space corresponding to the eigen- value — 2A for A = A^. Then we observe 
that along such a trajectory r _1 |V/| — > 2A. 

Let Pi < . . . < f3 n be the ordered set of eigen-values counted with multi- 
plicities. We also check that for e positive and small enough, the function 
|V/|i/-i(_ e ) has local extrema cj.(— s) < ... < c q (—e), for q > 1, which are 
continuous subanalytic functions of e. We also check that each such function 
c{-e)e~ 1 / 2 tends to 2/3 as e tends to for /3 = for some i. 
In other words the eigen-values 2Aj and 2fj,j are strongly connected with the 
limits at O of the local extrema of the restriction of |V/| to the (connected 
components of the) levels of the function /. 

5. Main result 

We use the notation introduced in Section [3l 

Non-oscillating hypothesis: We assume that \j\ does not oscillate at its 
end point O. 

We recall that plane gradient trajectories do not oscillate at their limit 
point, and that in higher dimension non-oscillating gradient trajectories exist 

as weii pa nam E3 . 

We begin with an elementary but key result about the monotonicity of 
functions along a non oscillating trajectory. 

Lemma 5.1. Let ifi : U —> M be a continuous semi-analytic function, an- 
alytic and non zero outside crit(/). The function has a finite limit at 
O, and up to shrinking U, the function il)^ is either constant or strictly 
monotonic on n U . 
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Proof. Assume is not constant. 

Assume there exists a sequence (sj)j converging to /(|7|), the length of 
I7I such that ifi' vanishes and changes sign at each Sj. Since 

f(s) = <V^,i//> 07(a), 

the subset | — y | n {(V^, Vf) = 0} has thus infinitely many connected compo- 
nents. Since clos({(V^,^j) = 0}) is a closed semi-analytic set and contains 
O, thus this fact contradicts the non-oscillation of |7[. □ 

The function s — > r = |7(s)| is strictly decreasing to as s —?■ l(T), 
following Lemma 15. 11 Moreover we have 

IV/I 

r(s) = o j(s) and r(s) — > — 1 as s — > ^(r). 

C'rJ 

So we reparameterize the integral curve 7 with r: 

l'{r) = ^;(7(0) = P^y u f) 70) = C 7W, 

where £ is the semi-analytic vector field appearing in the above expression. 
Another simple but key result is the following 

Lemma 5.2. Let a, b be continuous semi- analytic functions with {b = 0} C 
{d r f = 0}. Let v be the function defined as v(r) := (a/b) o^(r). Assume 
that v is bounded. Thus rv'(r) — > as r — > 0. 

Proof. First note that by the nature of 7(7") we will still have v' = a/b 2 
having the same form as v. From the proof of Lemma f5.ll we know that close 
enough to 0, the function v is strictly monotonic or constant. If v 1 —> 00, 
then rv 1 is still strictly monotonic close enough to 0. So if rv' does not tend 
to 0, integrating v' would contradict the fact that v is bounded. □ 

A first interesting consequence of Lemmas 15.11 and 15.21 

Corollary 5.3. lim^o Vf 07(f) = —v. 

Proof. First we write 7(7") = r{v + x( r )); where x( r ) = (pr — v) ~f(r). 
Since the vector field d r is semi-analytic, applying 15.11 we deduce that such 
that |x( r )| - > as r — > 0. Thus we deduce j'(r) = s'(r)u + j3{r) where 
P{ r ) = x( r ) + r x'( r )- From Lemma [5. II and Lemma 15.21 we get j3{r) = o(l). 
And to conclude we recall that s'(r) —1 as r —> 1, □ 

We recall that H(f) is "the normalized Hessian matrix" of /, and so the 
mapping x — > H(f)(x) continuous and semi-analytic. In particular thanks 
to Lemma 15.11 we deduce the following 

Lemma 5.4. There exists Ti € Sym(IR n )* such that 

tim\~f\3x^o H (f)( x ) = K- 



Note that this result can also be deduced from the results of [10] . 
This fact satisfied, the main result of this note is: 

Theorem 5.5. 

Hm HessOKx) • ujx) _ ^ H(f)(x)-u f (x) _ ^ 
| 7 |3x->o |Hess(/)0) • Vf{x)\ |7|3x-+o|fl'(/)(a:)-i//(aj)| 
Thus the vector u is an eigen-vector of the matrix H. 
The proof will consist of several elementary results. 

Proof. In the proof of Corollary 15.31 we had written 7(7-) = r(y + x( r ))> so 
that 7' = s'{r)v + (3 with f3 = x + r x'- Since 

x = d r -i> and x'(r) = -^V'/°7(?"), 

TO r J 

we can write /3'(r) = (a/b) o j(r) for a semi-analytic function v = a/b with 
i> _1 (0) = <9 r / _1 (0). Applying again Lemma I5T21 we deduce 

(5.1) 7 "(r) = oir- 1 ). 

Note that 

I V/l 

7"( r ) = v ' { r ) u f{l{ r )) + f 2 (^)(d^/ • z//) o 7(7") with i> (r) = o 7(f). 

We recall that v(r) — > —1 as r — > and so from Lemma |5.2| rv'{r) — > as 
r — >■ 0. We deduce from (15.11) 



(5.2) (di// • i//) o 7 (r) = o(r" 1 ). 

Since 

di// • u f = -—j- [Hess/ • v f - (Hess/ • v f ,v f )v f ] 



Equation (|5.2p now reads 

(5.3) r[Hess/ • £// - (Hess/ • z//, i>/)i>/] o j{r) = o(| V/|). 

The classical Bochnak-Lojasiewicz inequality applied to / and then to |V/| 
nearby the origin O gives 

\x\ ■ I Hess/ • V/| > C|V/| 2 for a positive constant C. 

or equivalently 

|x| • |Hess/ • Vf \ > C|V/| for a positive constant C. 
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Since r|Hess/| • Uf > C|V/| the only possibility for Equation (|5.3p to hold is 
that along [7I we get 

(5.4) Hess/ • Uf = (Hess/ • uj, Vf)vf + o(|(Hess/ • uj, Vf)\). 

Dividing both sides by |Hess(/)i^| and using Corollary 15 . 31 provides that the 
limit is either v or —v. Along [7I, the estimates / = —ar m d r + o(r m ) from 
Proposition 13.31 will give that it is —v. □ 

6. A FEW CONSEQUENCES 

For a given non-oscillating trajectory [ — ^ | with w-limit point O, we deduce 
from Proposition 13.31 Lemma 15. II and Lemma 15.21 that along (7! we have 

(6.1) V/ = -mar™" 1 ^ + o^" 1 " 1 ). 

The first un-expected result is the following one. 

Proposition 6.1. Let [7! be a non oscillating gradient trajectory of the 

function f with uj -limit point the origin O. Let m and a be as above. 

The real number m-a is an asymptotic critical value of the function r 1_m |V/|. 

Proof. Let us denotes by R the function r — > r 1_m |V/(7(r))|. Since R — > 
ma as r — > 0, Lemma 15.21 implies rR' — > as r — > 0. 
Since 

and — s> — 1 we get 

(rHess/ ■ u f - {m - l)\Vf\d r ,u f ) = o{r m ~ l ). 

Along I7I we know that V/ — ar m ~ 1 d r and so by Bochnak-Lojasiewicz 
r|Hess/ • Vf \ > Cr 171 " 1 . From Theorem 15.51 we also have |(Hess/ • vj, vj)\ ~ 
I Hess/ • Vf\, and so we necessarily deduce that 

rHess(/) -v f = {m- l)\Vf\d r + o(r m - r ). 

So we just deduce that along I7I the gradient vector field V|V/| is asymp- 
totically radial. Thus the result is proved. □ 

There are two consequences summarized in the next statements. The first 
one is already in the proof of Proposition 16.11 

Corollary 6.2. Along we find 

Hess(/) • Vf = m(m — l)ar m ~ 2 Vf + o(r m ~ 2 ) 

= —m(m — l)ar m ~ 2 d T + o(r m_2 ) 

= m(m — l)ar m ~ 2 v + o(r m ~ 2 ). 
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Let g be a homogeneous polynomial of degree p > 2. We recall that 

(6.2) p(p - l)g = (Ress(g) ■ rd r ,rd r ), 

For Corollary 16 . 21 and Lemmas 15. II and 15.21 we see that along I7I, we obtain 
an estimate similar to Equation ()6.2|) : 

m(m — l)f = (Hess/ • rd r , rd r ) + o(r m ). 

The second consequence is that the asymptotic critical value —a may be 
interpreted as an eigen-value under some additional hypotheses: 

Corollary 6.3. Assume that along \y\, r 2 h(j(r)) ~ r m ~ |/|. Then along 
I7I the matrix m ( m _\ym--2 Hess(/) has a limit at O in Sym(R n )* and the 
eigen-value in the oriented direction R + z/ is exactly —a. 

7. At infinity 

Even though the results of [5] were obtained in the semi-algebraic context, 
they still hold true for C 2 globally subanalytic function. 
Let / : R n — > R be a globally subanalytic real analytic function. We consider 
again the gradient differential equation: 

7(s) = f/(7(s)) with 7(0) = x and f(x ) < 0, 

But we are now interested in trajectories that leave any compact subset, 
that is, for s — > j(s) a parameterized solution, given any compact subset 
K of R n , there exists sk, such that 'j(s) ^ K for s > sk, or equivalently, 
|7(s)| — > +00 as s — >• 00. 

Assume that (7! is a non-bounded gradient trajectory, then | — y | leaves any 
compact. We still suppose that / is increasing along |7|. We may assume 
that the lim xe ^ x \^ +00 f(x) = 0. 

The main result of [S] is that the gradient conjecture is also true at infinity: 
Theorem 7.1 ([5]). There exists v € S™ -1 such that lim - — - = v. 

xe|7|,|a;|->+oo \x\ 

To prove this result, some estimates along the trajectory were necessary. 
We recall one that is relevant for this section 

Proposition 7.2 ([5]). There exists a positive rational number m and a 
positive real number a such that along \j\ 

f = — ar _m + o(r~ m ) and d r f = mar' 1 " 771 + o(r _1_m ) 

We need the following definition 

Definition 7.3. The trajectory \j\ is non-oscillating at infinity if for any 
globally subanalytic non-bounded analytic hypersurface H , there exists a ra- 
dius R such that \y\C\HC\ (R n \ B(0, R)) has finitely many connected com- 
ponents. 
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The main result in this context, using the notations introduced in the 
previous sections, is 

Theorem 7.4. Assume that \^\ is non-oscillating at infinity. Then 

1) lim xeWx ^ +00 H(f)(x) =H£ Sym(M n r 

2) The vector v is an eigen-vector of the matrix H. 

The proof will work along the same lines with a few modifications that 
we will describe. The first one is that there exists a Bochnak-Lojasiewicz 
inequality at infinity: 

Proposition 7.5 ([3]). If \x\ > 1 and \ f(x)\ < 1 then 

M • |V/| > C|/|, 

for a positive constant C. 

Then using it for the function |V/| we find again that for \x\ <C 1 and 
|V/| <C 1, there is a positive constant C. 

(7.1) \x\ ■ |Hess/ • V/| > C7|V/| 2 . 

Next we need the analogues of Lemmas 15.11 and 15. 21 

Lemma 7.6. Let g be a globally continuous subanalytic function outside of 
a compact subset ofW 1 . Assume that \^\ is non oscillating at infinity. Let 
v := g o As a germ at infinity, if v is not constant then it is strictly 
monotonic. 

Proof. If not it is contradiction to the non-oscillating at infinity of |'y|. □ 

Since the function 1/r is semi-algebraic, we can then reparameterize the 
solution 7 by r = |-y(s) | . 

Lemma 7.7. Let a,b be globally continuous subanalytic function outside of 
a compact subset ofM n , with {b = 0} C {d r f = 0}. Assume that [7) is non 
oscillating at infinity. Let v(r) := {a/b)o~f(r) and assume that v{r) — > c € M 
as r +00. Then the function r \— > r~ 1 v'(r) tends to as r — > +00. 

Proof. Similar to that of Lemmas 15.21 □ 

From these two Lemmas, and adapting the proof Corollary 15.31 it is easy 
to prove the following 

Corollary 7.8. Assuming that I7I is non oscillating at infinity, we get 

lim r ^ Vf l{r) = v- 

The proof of Theorem 17.41 goes along the same lines as those of the proof 
of Theorem 17.41 using Lemmas 17.61 and 17.71 Corollary 17.81 in a similar way. 

Along a non-oscillating trajectory I7I going to infinity, we find estimates 
similar to Equation (|7.2p . Proposition 16.11 Corollary 16.21 and Corollary 16.31 
namely 
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(7.2) V/ = mar- m - 1 d r + o(r~ m - 1 ). 

We also find 

Proposition 7.9. Let \j\ be a gradient trajectory of the function f non- 
oscillating at infinity. Let m and a be as above. 

The real number m-a is an asymptotic critical value of the function r m |V/|. 
Proof. It works as in the proof of Proposition 16.11 □ 
We eventually deduce the following estimates 

Hess(/) • Vf = —m(m — l)&r~ m ~ 2 Vf + o{r~ m ~ 2 ) 
= —m(m — l)ar~ m ~ 2 d r + o(r~ m ~ 2 ) 
= —m(m — l)a.r~ m ~ 2 v + o(r~ m ~ 2 ) 
m(m — 1)/ = (Hess/ • rd r , rd T ) + o(r~ rn ). 



8. RlEMANNIAN CASE 

Before getting in what really interest us, we will recall a few basic formulae 
of Differential Geometry. 

Let (M, g) be a smooth Riemannian manifold of dimension n. Let TM 
be its tangent bundle and 3£{M) be the C°°(M)-module of the smooth 
vector fields on M . Let (, } g be the scalar product coming with g and | — | g 
the associated norm. Let D be the canonical connection associated with g. 
Given X, Y, Z € J?T(M), the connection D satisfies the following equations: 

(8.1) [X, Y] = D X Y - D Y X = dY ■ X - dX ■ Y, 

(8.2) X ■ (Y, Z) g = (D X Y, Z) g + (Y, D X Z) S , 

(8.3) D X Y = dY -X + T(X,Y), 

Where T : 3£[M) x JT(M) -> 5£{M) is smooth, C°°-bilinear and symmet- 
ric. In local coordinates V is defined as usual with the mean of Christoffels 
symbols. 

Given a smooth function u, we recall that the gradient vector field V g ii 
is defined as the "dual" of Du for the pairing (,) g , that is for any smooth 
vector field X on M, 

(8.4) (Du)(X) = D x u = X -u= (V s u,X) s . 

The notion of Hessian of the function u as d 2 u is no longer accurate in 
this Riemannian setting. Since the connection D extends to tensors, the 
corresponding notion of Hessian in this setting comes from the following 
tensor 

D 2 u = D{Du) = (Ddu) : JT(M) x 3£{M) -> C°°(M) 

defined as 
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(Ddu)(X,Y) = (D x du)(Y). 

It is smooth, C°°-bilinear and symmetric as well. Using the formulae above, 
we also check that 

(8.5) (Ddu)(X, Y) = (D x V g u, Y) g = (D Y V g u, A) g , 
and deduce that 

(8.6) Vg|V g u|g = 2D Vg uV g u. 

We also look at Ddu as a C°°(M)-linear map X (M) JT(M), defined, 
with an obvious abuse of notation, as (Ddu)(X) = D x ^ s u. At each point 
x of M the tensor Ddu induces a bilinear symmetric form Hess g (u)(x) := 
(Ddu) x = T X M x T X M -)• R. We will call it the Hessian of u at x and 
obviously x — > Hess s (u)(x) is smooth. And we will also look at it as a linear 
endomorphism of T X M. 

Now we come to our topic. 

Assume that M = U is an open subset of M. n containing the origin O. 
Assume that the Riemannian metric g is analytic. 

Let / : U — > M be a real analytic function such that O € crit(/) the 
critical locus of /. We also assume f(0) = 0. 

Let us denote the unitary gradient vector field defined on U \ crit(/). 
Let us consider the following gradient differential equation: 

7(5) = J/f (7(s)) with 7(0) = xq and /(xq) < 0. 

We have picked some orthonormal coordinates at O so that g(O) is the 
Euclidean metric. 

If r is the Euclidean distance function to the origin O and r g the distance 
to the origin for the metric g, then, and r s /r — > 1 as x — > O. 

Remark 8.1. In a neighborhood of O, Bochnak-Lojasiewicz inequalities 
holds, namely 

rS|V g /|g > C|/| and so r|V g /| g > C'\f\, 
for some positive constant C, C . 

Moreover the results of [8] quoted in Section [3] still hold true, sticking 
with r in the estimates. 

Bochnak-Lojasiewicz inequality for the function |V g /|, then reads in a 
neighborhood of O 

r\D u sX7 g f\ g > C|V g /|g for some positive constant C. 

The operator Hess g (/) is now analytic and in the present coordinates, when 
we look at it as linear mapping, this gives us (symmetric) analytic entries 
hi j. We recall that the semi-analytic function h is h = (Y^ij hfj) 1 ^ 2 an d we 
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define H s (f) = h 1 Hess g (/), so that at each point x £ crit(/) we deduce 
H s (f)(x) is non zero. 

Now we can state the main result of this section in this general case. 

Theorem 8.2. Assume that I7I is a non-oscillating trajectory of the gradient 
field associated with f respectively to the analytic Riemannian metric g ; with 
w (| 7 [) = {O}. Then 

1) Along | 7 |, lhn xelMxH+00 HZ(f)(x) = H s G Sym(R^) \ 

2) 

Um _ lim H*(f)(x) ■ v f (x) _ y 

\-r\3x^o |Hess g (/)(x) • Vf{x)\ \~i\3x->o \H%(f)(x) ■ Vf(x)\ 

The vector v is an eigen-vector of the matrix T~L S . 

Proof. The first point is elementary. The proof of the second point goes 
exactly along the same steps as in the Euclidean case dealt with in Section 
[5j We nevertheless has to check the form of 7 in order to recycle the previous 
proof. 

j(s) = du^ ■ v^. 

Thus we find 

duf ■ uf = ^-[D^V s f - I>f , V g /) - (D uJ V s f, uf) g uf]. 

and we observe that r(z^, V g /) = 0(|V g /|). Thus applying Bochnak- 
Lojasiewicz for D u sV s f will allow to conclude again. □ 

Remark 8.3. 1) It is possible to find estimates as in Section^ but it is not 
clear what information some of them will really carry. 

2) In this context we will not speak of what is happening at infinity since the 
paper [5 J only dealt with the Euclidean metric. Moreover we would need to 
know some facts about the behavior of the metric g at infinity and it is clear 
what should be required to make things work. 

9. Remarks and comments 

First. What we have proved is not so much attached to a function than 
to a foliation, namely the foliation by the (germ at O of the) connected 
components of the levels of the given function. In particular, there is a 
unique "matrix" in the space of directions Sym(n) which is attached, at the 
limit point, to a given non-oscillating gradient trajectory. What this object 
encodes is still very unclear. We have clearly suggested that it is linked with 
the normal form of the divided gradient vector field, but also with the lines 
of ridges and valleys of the given foliation. 
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At this stage, it is not guaranteed that what we just suggested carries the 
relevant phenomena(on) to understand better what we have uncovered here. 

Second. Assume we are given an o-minimal structure A4 expanding the 
real numbers in which the semi-analytic sets are definable. We check that 
for C 2 plane function definable in A4 in a general C 2 definable Riemannian 
setting our main results is still true. In the Euclidean setting and at infinity 
we need globally definable to have it correct. 

If the o-minimal structure A4 is polynomially bounded, everything we proved 
in this note is still true for any C 2 definable function (globally definable 
at infinity), even the estimates The main reason for this is that Thorn's 
Gradient conjecture is still true in these categories [9]. Note that in the 
plane, gradient trajectories are still definable in an o-minimal structure: 
The Pfaffian closure of the initial o-minimal structure M. 

Third. Proposition 14.21 suggests that the limit % of our main result is linked 
with the reduction of the divided gradient vector field. Corollary I6.3l suggests 
that at such a point v the divided gradient vector field may be simple, 
but we have no idea whether it is true or not and then how to prove/dis- 
prove this. This is to be compared with the main result of [T], where is 
proved a resolution/reduction result at the end point of a subanalytically 
non-oscillating trajectory of real analytic vector field in a 3-manifold. 
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